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Abstract 

A calculation of the hadronic matrix elements for K — > tttt including tt-tt 
rescattering effects in a dispersion integral is presented. I study the dependence 
of the results on the matching scale u and use them to calculate the CP-parameter 
(ef/e). I find improved stability on the matching scale and good agreement with 
the experimental results. 



During the past few years there have been new developments on the problem of direct 
CP-violation. On the theoretical side there are predictions which agree with the 

new experimental results || [7] and some calculations which are still far away from the 
data. The agreement between theory and experiments must be considered a success of 
QCD since the dominant contributions to {si/e) come from the effective Hamiltonian gen- 
erated through renormalization of the weak interaction. The renormalization coefficients 
were computed by two groups || || and agree with each other. 

It is clear now that the original calculations of the matrix elements, done in the Vacuum 
Saturation Approximation (VSA), must be supplemented by rescattering corrections of 
the two pions, which increase (Qe)o an d decrease (Qs)2- The rescattering corrections 
were described in an earlier calculation |2j which made the ratio (e//e) positive. A new 
and improved calculation in the chiral theory of pseudoscalar mesons 0, || and the chiral 
quark model obtained results consistent with large values for (st/e). In this article I 
introduce a new method for calculating the rescattering corrections, which may also be 
useful for other low-energy calculations. The results indicate that tt-tt rescattering brings 
important corrections and allow to compute the imaginary parts of the amplitudes. 

The operators occurring in the effective Hamiltonian are local operators which satisfy 
analyticity and threshold conditions of field theory. Consequently their matrix elements 
can be written as integrals over their singularities. There are two types of contributions 
to each matrix element: tree diagrams and loops. The four dimensional integrals for the 
loops can be rewritten as one-dimensional integrals over the energy flowing through the 
K-meson. We split the integrals in two regions - from 4m 2 to /i 2 and from /z 2 to infinity. 
The higher energy region is computed in QCD and its discontinuity is related to the 
anomalous dimension. The low energy part can be represented by an effective theory or 
by experimental data. The final formula for the amplitude at a low energy point o = m 2 K 
can be written as 
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Here a is the lowest order contribution, which can be substituted by the vacuum saturation 
or another improved approximation. In the last equation we factorized two terms and 
identified one of them with the series generated by QCD. We emphasize that we do not 
need a subtracted dispersion integral because for the high-energy region we use directly 
the QCD realization of the theory. This is a useful expression for the low-energy value 
of the amplitude and may have applications beyond the cases described in this article. 
In addition, eq. (1) provides a matching between the two regimes at the cutoff scale fj,, 
which can be varied from 0.8 GeV all the way up to 1.8 GeV where the validity of QCD 
is more reliable. 

In the K — > 2-k decays we have many matrix elements of the operators, generated 
through QCD, but very few observables. To improve the situation we shall try to calculate 
the amplitudes and their strong phases. For the A 2 amplitude we use the value 

A(K + -> tt+tt ) = (1.837 ± 0.002) • 10~ 8 GeV (2) 

and its phase is given, according to Watson's theorem, by the phase extracted from n-n 
scattering. In fact we will calculate the phase for each matrix element and compare it with 
the 7r-7r phase shift. The imaginary part for each rescattering is given through unitarity 
by the cut diagram shown in figure 1; where the square is a weak vertex and the circle a 
strong vertex. We obtain 

Im(nn, I = 0\Q 6 \K°) = (ttvt, / = 0\Q 6 \K°) VSA \ ^"Jf* (l " ^) ^ (3) 
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I shall use these equations for computing the imaginary parts throughout this paper. Using 
them we computed the numerical values for Im(Qs)2 in Table (1) of reference The 
imaginary part of (Qe)o was not included in [Q because it is higher order in the expansion 
described in that article. The form of the equations indicates that the imaginary part 
consists of two multiplicative factors: a weak vertex and the n-n scattering amplitude 
characterized, in our case, by two isospin states I = 0, 2. 

We proceed now to calculate the real part of the amplitudes. For this calculation we 
need the imaginary part for values of the center-of-mass energy squared, s, in the range 
4m 2 < s < yU 2 . The proof of Watson's theorem holds for the matrix element of each 
operator and their imaginary parts are given as 

lmA l = Jm(vr7r, I\Qi\K°) = ((tttt, I\Qi\K°}\ ■ smSj^s) (5) 

with Sq the experimental phase shifts for isospin / = 0, 2 pion-pion scattering. The 
imaginary part is given by this formula in the elastic region and I adopt this form beyond 
the elastic region using the experimental phase shifts. For the magnitude of the matrix 
element we can take the low energy contribution, mentioned earlier, but we are free to 
introduce a weak energy dependence; for instance the variation introduced by the real 
part obtained through the dispersion relation. Let us denote by 

A Q ,,(s) = (tttt, I\Qi\K°). (6) 

Then 

ReA Qi>I (a) = \a Qi + - f ((tttt, I\Qi\K°)\ da) (7) 

We shall assume that the absolute value of the matrix element is a slowly varying function 
of energy over the region of integration and we use the experimental values for siuSq to 
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perform the principal value integral. We define the functions 
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and calculated the values presented in table 1. These terms bring in a correction to the 
tree level contribution. For the integrations we use the experimental phase shifts from 
references ||10||-||13||. Data for 5q exist up to the energy of 1.5 GeV and we performed the 
integral up to this value. For 5® the data show clearly the p-resonance. They extend to 
1.8 GeV and I give the additional values in the table. The functional change of the phase 
shift given in flTT] produces an fo(fi,mk) which is almost constant for /i > 900 MeV. The 
extended ranges presented in table 2 allow us to study the dependence of ReA Qi I on the 
matching scale \i. As an additional test, I introduced a 20% energy dependence on the 
magnitude of the matrix element and carried out the integration. The results differ only 
by one or two units in the second decimal. 



As a first test of the approach we consider the A 2 amplitude. Its Born contribution is 
known 
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^ ^ V ud V: s {zx + z 2 )(Q t } 2 = 2.74 x 10- 8 GeV. (9) 



fj, in GeV 


f 2 (fi,m k ) 


fo(fi,m k ) 


0.7 


-0.09 


0.34 


0.8 


-0.12 


0.50 


0.9 


-0.17 


0.60 


1.1 


-0.23 


0.71 


1.3 


-0.29 


0.64 


1.5 


-0.34 


0.55 


1.7 




0.55 


1.9 




0.60 



Table 1: Numerical results for the principal value integrals 
Including next the unitarity corrections, at /i = 0.9 GeV we obtain 
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2.74(1- 0.17 -i0.20) • 10" 
(2.27-i0.55) ■ 10" 8 GeV 
-0.24, 6 = -13.7°. 



(10) 
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We note that the magnitude of the calculated amplitude is larger than the measured value 
by 30%. The phase of the amplitude has the correct sign and it is slightly larger than 
the experimental value. The experimental values reported in the articles [|H|~[|HI var Y 
among themselves. An approximate value is —10.5 ± 1.6°. A more accurate value was 
obtained through a dispersion calculation for n-n scattering JTIJ, but the error quoted is 
very small and is perhaps an underestimate. Finally, the dependence of the amplitude on 
the matching scale \i is small, as discussed for the other two matrix elements below. 

The calculation for Aq g 2 proceeds along similar lines. We obtain the imaginary part 
from eq. (4) and the real part from eq. (7) and table 1 



A Qa2 (m k ) = (vrvr, 1 = 2|Q 8 | j FT°)vsa(1 - 0.20-^0.20) yM 



at fj, = 1.0 GeV and with §, § 
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have the numerical values 



The renormalized couplings L 5 and L 8 are defined in references 



(12) 



(13) 
and 



L 5 = 2.07-10" 



and L 8 = 1.09 • 10~ 3 . 



Substituting the numerical values for m s (1 GeV) = 150 MeV we obtain 



^38,2 (m k ) = (0.37-i0.09)j/ 8 (/i) GeV 3 . 



(14) 



(15) 



The phases in eqs. (10) and (15) are the same, which is an attractive property of the 
dispersion relation, i.e. all I = 2 matrix elements have the same phase. A similar property 



holds for the 1 = matrix elements. Finally, we can investigate the dependence on the 
matching scale \x which appears on the factor r 2 , through the running mass of the strange 
quark, and the Wilson coefficient y^i)- In a specific regularization scheme the product 
r 2 ys(fi) is stable between 1 and 2 GeV, varying by less than 5%. The variation among 
the three regularization schemes in the same energy region flTJj is at most 15%. 



The unitarity corrections to the matrix element Qq are controlled by the phase shift 5q 
which is positive, giving a positive correction for the real part of the matrix element. The 
increase of the matrix element is given by the function /o(/f) in table 1, which is stable 
above 900 MeV. At n = 1 GeV 

Aq 60 = (HQ 6 |ir >vsA(l + 0.60 + «0.40)?/ 6 (/x) 



with § 1 



(-0.56 - zO.14) y 6 (ji) GeV J (16) 



MQ 6 |^ >vsa = -4v/3r 2 L 5 ™ l ™ l . (17) 

-T-7T 



Again the product r 2 ?/6(/-0 is very stable for 1 < // < 2 GeV; the variations mentioned 
in the previous paragraph for r 2 y 8 (fi) again hold. The phase of this amplitude is positive 
and equal to +14° which is approximately half the experiment phase shift at yfs = m k . 

With the values derived already we can compute the parameter (ef/e). The standard 
derivation leads to the expression 
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with 

tt 2 = |£ i2/i Gu)<Q,} 2 | (20) 

Y = 7T --7T 2 (21) 
0J 

and Q m , ~ 0.25 ± 0.05 being the isospin breaking in the quark masses (m u ^ rrid). 
The absolute values originate from the fact that the phases of strong origin were already 
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extracted in the calculation of ei/s. The overall factor is precisely known ^ \ £ \r £ a 
GeV -3 and the Cabbibo-Kobayashi-Maskawa factor was recently estimated |i"4fl . 



ImX t = (1.38 ±0.33) x 10" 



346 



(22) 



This new value is a large improvement over the values reported in the early 90 's and leads 
to a large reduction of the uncertainties. A second reduction of uncertainties comes from 
the weak dependence of the amplitudes in eqs. (12), (15) and (16) on the matching scale. 
We summarize in table 2 the uncertainties for (s//e) originating from two sources. 



source 


A(e,/e) 


Im Xt 
Matching 


±25% 
± 20% 



Table 2: Theoretical Uncertainties 



Another uncertainty comes from the strange quark mass which enters the calculation of 



the matrix elements through the factor r 



_2m| 

m s +fh 



-. Values for the running strange quark 



mass have been computed by various methods and vary considerably. Older estimates 
IS] and QCD sum rules [TI| give higher values m s (l GeV) ~ 150 ± 55 MeV at the usual 



renormalization point fi = 1 GeV; while recent values from lattice calculations |Ty| give 
smaller values ~ 110 ± 20 MeV at /x = 2 GeV. It is customary to adopt the range [[[ 



m s (lGeV) = 150 ±25 MeV. 



(23) 



For central values of the parameters at 1 GeV and keeping the errors from table 2 in 
quadrature, I obtain 



ef/e = (4.2 ± 1.3) ■ (0.0320) ■ 10" 



'15.0 ±4.8) ■ 10" 



(24) 



where the number 0.0320 comes from (Qq)o including Q vv , minus the contribution from 
(Q%l2- I kept only these two operators and found out that the contribution from the 
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1 = 2 amplitude is only 20% of the (Qq){1 — Vt m ,) term. This value for the ratio is in 
good agreement with the average experimental value {si /e) = (21.2 ± 4.6) ■ 10 -4 . 

I have shown in this article that the VSA for the matrix elements together with 
rescattering corrections lead to values of (e//e) which are consistent with the experimental 
measurements. The unitarity corrections improve the stability on the matching scale [i. 
The specific values of the phase-shifts increase (Qe)o and decrease (Qs)2 making the 
difference in the function Y positive definite [Q - a feature which is maintained in many 
calculations f2|-||, |L9| . Two very recent articles [|l^, [HJ use dispersion relations for the 
K-meson decay amplitudes. They differ in several basic respects from the present article 
and the interested reader can study them for comparison. Finally, it will be interesting 
to test if lattice calculations [p0|, [2l[] , which fullfil unitarity corrections, still give different 
results. 

A more extensive exposition of this work including numerical studies and other appli- 
cations will be presented in the future. In particular, I wish to study the various contri- 
butions to the A/ = \ rule and compare them with the results of previous calculations 
El-El. 
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